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I. INTRODUCTION 



In the years 1893-96, the Norwegian explorer Fridtjof Nansen, while he were traveling 
in the Arctic region, noticed the drift of surface ice across the polar sea making an angle 
of 20 to 40 degrees to the right related to the wind direction. Nansen has advanced an 
explanation, that in addition to the wind's force, it was also actuating the Coriolis force. In 
1905, Vagn Walfrid EkmarP introduced a theory of wind currents in open seas showing that 
sea current change of direction with depth due to Coriolis force appearing in a rotating based 
coordinate system. Besides attracting considerable interest in geophysical flow problems^, 
these discoveries stimulated further investigation in such fields as magnetic geodynamics, 
binary stars, new-born planetary systems, trying to unveil the important problem of angular 
momentum transport, and as well, this work. 

The problem of enhanced angular moment transport in accretion diskP, and the break- 
down of Keplerian rotation, as well as the removing of angular momentum from a vortex 
due to moving spiral waves, playing an important role on the total angular momentum bal- 
ance of the core and the intensification of a tropical cyclone^, are all examples of problems 
that demand a clear understanding of the dynamics of electromagnetic-gravitational rotat- 
ing systems. Furthermore, special attention must be dedicated to the role of flux of angular 
momentum, and its conservation. Notwithstanding all attempts to take an accurate account 
of the transport angular momentumP, it lacks an additional equation of conservation (besides 
continuity, momentum and energy equations), which is related to local density of angular 
momentum, and its related flux. In order to develop a logically closed foundation, the equa- 
tion for the angular momentum balance must be included and, to our best knowledge, these 
problems were treated by Curtiss^ and Livingston and Curtiss^, and they are also addressed 
in the present work but using a variational method. 

In this article, we develop a standard technique for treating a physical system on the 
basis of an information-theoretic framework previously developecP^, which constitutes an 
entropy maximizing method leading to a set of two first-order differential equations revealing 
the same formal symplectic structure shared by classical mechanics and thermodynamics^^. 
Our method bears some resemblance with isoentropic but energy-nonconserving variational 
principle proposed by R. Jackiw et al. 13 allowing to study nonequilibrium evolution in the 
context of quantum field theory, but with various classical analogous such as the Schrodinger 
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equation giving rise to reflectionless transmission. 

The plan of the article is as follows. In Sec. II, we include the extended mathemati- 
cal formalism to nonequilibrium information theory. In Sec. Ill, we analysis the equilib- 
rium and stability of a rotating plasma. In Sec. IV, we apply our formalism to angular 
momentum transport, obtaining in particular the Umov-Poynting's theorem for rotating 
electromagnetic-gravitational systems (e.g., rotating plasmas, magnetic geodynamics, vor- 
tex motion and accretion disks in astrophysics), whose applications might contribute to 
clarify still poorly understood phenomena. 



II. MATHEMATICAL PROCEDURE 

Let us consider here a simple dynamical system constituted by a set of N discrete inter- 
acting mass-points (a = 1, 2, N) with xf and vf (i — 1, 2, 3; a — 1, N) denoting 
the coordinates and velocities of the mass point in a given inertial frame of reference. The 
inferior Latin index refers to the Cartesian components and the superior Greek index dis- 
tinguishes the different mass-points. 

The gravitational potential 0*°) associated to the mass-point a is given by 

>) = G r m(a) m 

^ U 2^ | x (a) _ x (/3)|' W 

0=1 1 1 

with G denoting the gravitational constant and x( Q ) and x^) representing the instantaneous 
positions of the mass-points (a) and (/?). Yl denotes summation over all the particles of the 
system. Energy, momentum and angular momentum conservation laws must be verified for 
a totally isolated system: 

TV 

E = Y,E ia) , (2) 



a=l 
N 



£p (a) > 0) 



a=l 



the particles total angular momentum (sum of the orbital angular momentum plus the 



intrinsic angular momentum or spin) 

N 

L = Y, L (a) = $^([r (Q) x p (Q) ] + J (Q) ). (4) 



N N 



a=l a=l 



In the above equations, is the position vector relatively to a fixed frame of reference 
K, p (Q) is the total momentum (particle + field) and L*^ is the total angular momentum 
of the particle, comprising a vector sum of the particle orbital angular momentum and the 
intrinsic angular momentum J (e.g., contributed by the electron spin and/or nuclear spin, 
since the electromagnetic momentum is already included in the preceding term through 
pO)) The maximum entropy principle introduces Lagrange multipliers from which, as we 
will see, ponderomotive forces are obtained. 

It is necessary to find the conditional extremum; they are set up not for the function S 
itself but rather for the changed function S. Following the mathematical procedure proposed 
in Ref.^ the total entropy of the system S is thus given by 

S = t {s^ - _ W _ q M V i*) + ^(aW . V W) _ u<A 

a=l * ^ ' 

JV 

+ (a-p (Q) + b- ([r (Q) xp( fl )] + j( a ))} =J2 &{a) - ( 5 ) 

a=l 

where a and b are (vectors) Lagrange multipliers. It can be shown that v re ; = aT and 
lj = hT (see also Ref.^. The conditional extremum points form the dynamical equations 
of motion of a general physical system (the equality holds whenever the physical system 
is in thermodynamic and mechanical equilibrium), defined by two first order differential 
equations: 



OS 

> ► canonical momentum; (6a) 



dp(°0 



dS 



1 1 <9v( a ) 

— — V r ( Q )£/^ — —m^ — - — > ► fundamental equation of dynamics. (6b) 
or^' T T at 

Eq. [6b] gives the fundamental equation of dynamics and has the form of a general local 
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balance equation having as source term the spatial gradient of entropy, V a S > 0, whilst 



Eq. 6a gives the canonical momentum (see also Eq. 12). At thermodynamic equilibrium the 
total entropy of the body has a maximum value, constrained through the supplementary 
conditions |2j [3j and |4j which ensues typically from the minimization techniques associated 
to the Lagrange multipliers. In the more general case of a non-equilibrium process, the 
entropic gradient must be positive in Eq. 6b, according to Vanderlinde's proposition ^1, 



a condition required for the gravitational force to exist. However, new physics may be 
brought by the set of two first order differential equations, related to the interplay between 
the tendency of energy to attain a minimum, whilst entropy seeks to maximize its value. 

In non-equilibrium processes the gradient of the total entropy in momentum space mul- 
tiplied by factor T is given by 

<9p(") [ m (°0 TO («) L J J W 



so that maximizing entropy change in Eq. 6a leads to the well-known total (canonical) 
momentum: 

p (a) = m (a) Ve + m (a) [ w x r («)] + A. (8) 

The above formulation bears some resemblance with Hamiltonian formulation of dynamics 
which expresses first-order constraints of the Hamiltonian H in a 2n dimensional phase 
space, p = —dH/dq and q = dH/dp, and can be solved along trajectories as quasistatic 
processes, revealing the same formal symplectic structure shared by classical mechanics and 
thermodynamics^^. 

The internal mechanical energy term, UmL, appearing in Eq. [5] may be defined by: 

N 

Uj£l = m (c V Q) (r) + ^ '■ ( 9 ) 

p=i 

Considering that by definition of thermodynamic temperature, dS^ / dU^ = 1/T^ a \ then 
it follows 

V r(ci) U (a) = -m^V<p ia) -^-Vv {a) -V r(a) (j^j -g (Q) VV (Q) +g (a) V(A Q -v^). (10) 
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Eq. 10 contains the particle's self-energy and the particle interaction energy for the gravita- 
tional and electromagnetic fields, but it may also include other terms representing different 
occurring phenomena (energy as a bookkeeping concept), such as terms included in Eq. [9j 
We may here recall that the entropic flux in space is a kind of generalized force X a H3EE3 
and it can be shown that the following equation holds: 

TV r(Q) 6 (a) = -q {a) V ria) V {a) + g (a) V r(Q) (A (a) ■ v (a) ) 

+ ro<«M«> ■ Vv<«> - V rW (y^?- - u ■ J (a) ) , (11) 

Now, it make sense to write the fundamental equation of thermodynamics under the form 
of a spacetime differential equation: 

rV5+ E m(a) ir = E vu<a) - ( 12 ) 

a a 

Taking into account the convective derivative, dv^/dt = dv^ /dt + v^ -Vv^, after some 
algebra we obtain: 

m (o) — - = -TV r(Q) 6 Q - m (Q) V0 (a) - g< a >VV (o) 
dt 

+ gWV(AW • vW) - V rW (^f - cj ■ 3^ + f2. (13) 



For conciseness, the term [7^ a ^ now includes all forms of energy inserted into the above 
Eq. [5j On the right-hand side (r.h.s.), the first term must be present whenever the mechan- 
ical and thermodynamical equilibrium conditions are not fulfilled; the second term is the 
gravitational force term; the third and fourth terms constitute the Lorentz force; the fifth 
term is a new term which represents the transport of angular momentum; while the last 
term represents other external forces not explicitly included and actuating over the particle 
(a). 
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A. Examples: 



1. Extended fundamental equation of dynamics 

The present formalism was applied in a previous worlpf, and therein we obtained the 
ponderomotive forces acting on a charged particle. For a neutral particle (body) in a gravi- 



tational field, Eq. 13 points to a kind of extended fundamental equation of dynamics for a 



given species (a) at equilibrium and at a given point of space-time (Eulerian description): 



dt 



-m^W a) ~ V r(a) - u • J (a) ) . (14) 



Eq. [14] has a new term because the body possess an intrinsic angular momentum. In a 
non-rotating frame of reference we must put u = 0, and from the work-energy theorem it 
is obtained the total mechanical energy of the system, £ mec = K + U + J 2 /2I C . This is the 
common approach in classical mechanics. We are interested in the effect of a given force 
in a given point of space-time, not in its effect along the particle trajectory. It is worth to 



point out that Eq. 14 was obtained through a variational procedure in contrast to the usual 



conservation theorem used, for example, in Ref. | 18 | 19 i 



Included in the internal energy term are the interpressure term, see Eq. [9] (we consider 
here a homogeneous and isotropic fluid). The above described framework (see also Ref.^ 
for additional information) leads us to the well-known hydrodynamic equation for a given 
species (a): 



0=1 



- v ^> -u-J^ + Te^). (is) 

Here, in the r.h.s. of the Eq. 15 we introduce explicitly external forces terms, eventually 
present in open systems. 

Using the following correspondence from particle to fluid description 

J2^ {a) -> J d 3 x' Pv {x'), (16) 



and, as well, a analogue relationship for the electric charge 



,(«) 



<i 3 x p(x'), 



(17) 



we can can rewrite Eq. 15 under the form of the Euler (governing) equation: 



dt 



Pv-JT = -/'rV r O - V r p - V r $j - V r /. 



(18) 



Here, as usual, the total interparticle pressure term (e.g., RefP^) is given by: 



N 



p(r) = ^m( Q )^0^(r). 



(19) 



/3=1 



To simplify, we introduce a functional integral in the form of an intrinsic angular momentum 
energy density (comprising the "interaction energy term", U) ■ 3), $j: 



j(ap 



-a/ - J (Q) - (AF) 



(a) 



(20) 



considering that the intrinsic angular momentum density refers to a given blob of fluid (with 
inertial momentum J, a measure of the local rotation, or spin, of the fluid element), and 



its associated free energy / = f — Ts (per unit of volume). Eq. 18 also means that 



function (the field integral of r^ a ') is constant along the integrals curves of the space 
field r' a l The gradient of the free energy / of the out of equilibrium state is the source of 
spontaneous change from an unstable state to a more stable state, while performing work. 
For example, a common source of free energy in a collisionless plasma is an electric currentpl 
in a magnetically confined plasma, several classes of free energy sources are available to 
drive instabilities, e.g., relaxation of a non-Maxwellian, nonisotropic velocity distribution^. 
At this stage, it is worthwhile to refer that our procedure includes the effect of angular 
momentum (through Eq. [5]), as it should be in a consistent theory, according to Curtiss^. 



Using, furthermore, the mathematical identity: 
V(A (Q) ■ v (a) ) = (A (a) • V)v (a) + (v (a) ■ V)A (Q) + [A (a) x [V x v (Q) ]] + [v (Q) x [V x A (a) ]], (21) 
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we obtain after some algebra the following expression: 

V (A (a) • v (Q) ) = [u x A (Q) ] + [v (Q) x B (a) ] . (22) 

Here, B = [V x A]. In addition, we may notice that the following equality holds: 

(A (a) • V)v (a) = [u x A (a) ], (23) 

where 

A M = ^g^ — , (24) 

/3=1 T °$ 

denotes the vector potential actuating on (a) particle due to all the other particles, and the 
vorticity is defined by 

= [V r(Q) x v^] = 2u/ a) . (25) 

Therefore, it follows the general equation of dynamics for a physical system (Lagrangian 
description) : 

ctv 

p— = pE + [ J x B] — V0 — Vp + p[A x u]. (26) 



The last term in the r.h.s. of the Eq. [26] is a new term that represents a kind of topological 
spin vector^, an artifact of nonequilibrium process. As it will be shown subsequently, 
the topological spin vector plays a role in plasma arcs, and as well in magnetocumulative 
generators^, and suggests a new method to obtain the helicity transport equation^. 



2. Rigid body rolling down an inclined plane with holonomic constraint 



We may check Eq. [14] in a standard example of classical mechanics: a rigid body of mass 
M rolling down an inclined plane making an angle 9 with the horizontal (see, e.g., p. 97 
ofRefP). Eq.p can be applied to solve the problem, with u = (there is no rotation of 



the frame of reference) and considering that only the gravitational force acts on the rolling 
body, with inertial moment relative to its own center of mass given by I c = f3MR 2 . Hence, 
we obtain: 

Mx = Mg sin 9 - d x w. (27) 



Here, w = 2/ ■ Assuming that the x-axis is directed along the inclined plane, and con- 
sidering that the angular momentum relative to the rigid body center of mass is given by 
J c = I c u', with oj' = d6/dt, and noticing that dx = v x dt (holonomic constraint), it is readily 
obtained 

duj' uj' 

M'x = Ma x = Mg sin 9 - I C J — - = Mg sin 9 - (3 MR 2 a. (28) 

v x dt ' co'R 

Since a = a/R, then it results the well-known equation 



gsin9 



(29) 



III. EQUILIBRIUM AND STABILITY OF A ROTATING PLASMA 

Extremum conditions imposed on entropy or internal energy, does not only provide crite- 
ria for the evolution of the system, but, in addition, they determine the stability of thermo- 
dynamic systems, at equilibrium. Further, it has been showrr^ that a state of mechanical 
equilibrium can be reached, if the entropy increases with distance: 

TL r S > 0, (30) 

where L r is the Lie derivative along the vector field r acting on scalar S. Therefore, we can 
search another extremal condition through the general expression 



= - 9 ( °»V r ,^'»» 4- ,<«>V r ,.,(v«-> • AW) - m <«>^ - V rW (<g£ - « ■ J«>) . 

(31) 

According to Noether's theorem, the total canonical momentum is conserved in a closed 
system, and hence we can state the closure relation: 



N 



[P (a) + <? (Q) A (a) ] = 0. (32) 



a=l 



It can be shown (see RefP^ for details) that the relationship that prevails in equilibrium 
on a rotating plasma is given by 
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g («) E M + gW[v x B (Q) ] = m (Q) V r(Q) (a) + m (a) ^ V r(Q )0 (Q,/3) - g (a) [A (a) x w^)]. (33) 



We may develop the correspondence drawn from Eqs. 16- 17 and to obtain the general 
condition of equilibrium of a rotating plasma, in the presence of gravitational and electro- 
magnetic interactions (e.g., RefP^): 

pE + [ J x B] = V$ + Vp — p[A x u]. (34) 

Hereinabove, to simplify the algebra, we took the averaged angular frequency u, for the 
all system. Embedded into Eq. [34j it is present the vector potential A, on the same foot- 
ing as the E and .B-fields. The importance of the vector potential when compared to the 
others, depends on the kind of I?-field prevalent on the system; for instance, if the I?-field 
is homogeneous, the vector potential field predominates on the region near the axis, since 
puA/JB ~ 1/r 2 . The topological spin vector term is fundamental, since it produces work, 
responsible for the system angular momentum modification, producing a rocket-like rota- 
tion effect on the plasma. The theoretical framework delineated here, may help to clarify 
problems related to rotating-plasma systems^^, and controlled thermonuclear plasma con- 
finement 



A. Bennett pinch 

The compression of an electric current by a magnetic field, the z-pinch effect, can be 



studied on the basis of Eq. 34 , which gives the condition for dynamic equilibria. Let us here 
assume a typical geometry for an infinitely long axisymmetric cylindrical arc (Fig. [TJ, with 
axial current density J z = J z {r). Since the current density is assumed constant, Maxwell's 
equation in the steady state yields the azimuthal component Bg = /j,oJ z r/2, for r < R, with 
R the outer boundary of the cylindrical arc. The vector potential is purely radial, given by 



A z {r) = —fioR 2 J z /A, for r < R, and the Coriolis term plays no role. We can write Eq. 34 
under the form: 

- J A = % (35) 
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FIG. 1. Geometry and vectorial fields in the Bennett pinch generated by an axial current J z 
creating a toroidal field B^. If, instead, we consider a vacuum arc discharge with radial current J r 
and magnetic field B z , we may have rotating arc with angular velocity u. 



and it follows that 



p(r) 



R tr ir= W' <H2 - r2) 



(36) 



which is a well-known result. 



The interaction between vacuum arcs and transverse magnetic fields, is used in switching 
devices (see e.g. RefsJ^ESl). We can consider here instead, a coaxial configuration with a 
cathode on axis with a stabilizing magnetic induction field B, directed along the symmetry 
axis, and an arc current density J flowing radially (and assuming a "filamentary" current 
with radius R', A r = —fioR' 2 J r /A, with /j, the permeability of the vacuum). In this case, 



we may apply Eq. 34 to obtain the pressure differential, from the axis to the wall (at R): 



Ap(r = R) = 2-kR 



■J r B z + I p c I T[R 2 ^-J r uj 

4:71 



A9. 



(37) 
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. . Pq B z . . 

I p c \> — — retrograde rotation (38a) 
I Pc |< ~^r w > amperian rotation. (38b) 

Here, S denotes the filamentary current cross-section, and p = p c . For negative charge 
carriers p c = — \ p c \= — en e , and we obtain an amperian (clockwise) rotation for high 
magnetic fields, and relatively weak arc currents, or instead, a retrograde rotation for higher 
intensity arcs (higher S) and small transverse magnetic fields, in agreement with experimen- 
tal evidence (e.g., RefP^). We may see in this effect the compromise between two different 
tendencies, the energy which tends to a minimum and the entropy which tends to a maxi- 



mum. From Eq. 37 it is obtained an expression for the spot velocity in a transverse magnetic 
field: 

Air 1 f e i B Z R\ 

v = uR=—- r + -S-J- (39) 

p | p c | \2na c S J 

Here, we have made use of the Bernoulli relation, Ap = eoE 2 /2, and of the constitutive 



equation J r = <j c E, <j c denoting the plasma electrical conductivity. Although Eq. [39] is not 
self-consistent, it shows that the force term [J x B], is not playing alone its key role, and the 
spot velocity depends linearly of the arc current (see, e.g., Ref.^). A better understanding 
of this phenomena is crucial, since arc discharges are powerful generators of non-equilibrium 
atmospheric pressure plasma. 



IV. FLUX OF ANGULAR MOMENTUM 

We now address the transport of angular momentum, a phenomena of huge importance 
in several phenomena, such as in the working mechanism of an accretion disk, the formation 
of a tornado, or the atmospheric circulation. 

The equation of conservation of angular momentum can be obtained from the following 
equation: 
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+ g W V r(Q) (vW ■ AW) - V r <„> - u* ■ J<°)) 



(40) 



Multiplying the previous Eq. 40 by the particle velocity ?4 (z = x,y,z), and after rear- 
ranging the terms, we obtain: 

9^ J+ l 1 1 9 Z-^i Z-^i | x (a) _ x (/?)|3 



2^^te 1 J |x(«)-x(«| 3 2^- W 4 ^ j 

a p^a u 1 1 a 



1 \^ (q) 
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<9x 



(a) 



J(a) |2 

2J( Q ) 



(a; Q ■ J (a) ) 



,(«) 



<9x 



(a) 



(41) 



Eq. 41 can be written in a more comprehensive form, if we separate the terms from differ- 
ent contributions. For this purpose, let's define the total kinetic energy by means of the 
expression 

j = lj2m {a) \v {a) \ 2 , (42) 



and let us denote by 



9 ^ ^ | x (a) _ X (/J)| 



(43) 



the overall gravitational energy. Similarly, the total electrostatic energy is given by: 



u a fi^a 1 1 



(44) 



Following now the Umov-Poynting procedure, it can be shown that the previous Eq. 40 
reduces after some algebra to the following one: 

d '3) = - ip„ + \ J2*> [a) x t A(a) x W H 
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re* 1 2 



2/( Q ) 



- (w • j (q) ; 



,(«) 



9©(a) 



(45) 



As we did before, it is convenient to introduce a new physical quantity, representing the 
rotational energy: 



q3 



rot 



(46) 



Then, the above Eq. 45 , can be written under the more general form: 



d 
di 



(3 + Vp g + Vp v +Vp 



rot j 



-V- V t 



21 



— (u ■ 3) — AF 



(47) 



Here, v represents now the velocity of an "element of fluid" , and the last term in brackets 
is an average local value. A simple analysis of the Eqj47] shows that the system is stable in 
rotatory motion, provided that the following condition is satisfied: 



21 



(w-J) 



AF 



< 0. 



(48) 



The last term is the free energy per unit of volume, F = Fq — TS. Remark that Eq. [48] is 
consistent with Gibbs distribution in a rotating body (see, e.g., Ref. ™, which means that 
the radial flux of energy must be positive, flowing out radially from the system's boundary. 
In addition, we notice that the equilibrium of a gravito-electromagnetic system depends 
upon its mechanical rotational properties, but also on the free energy available, linking 
intrinsically any mechanical process to thermodynamic variables, and opening options for 
possible unconventional mechanisms to control instabilities. In the domain of astrophysical 
plasmas, gravity and rotation usually dominate over the magnetic field of force, being crucial 
for the development of instabilities. 



From Eq. 47 we may see that equilibrium ensues (neglecting thermal and configurational 
effects) when the velocity of rotation of the fluid satisfies the local condition (e.g., RefP^): 



dfl(r) 
dr 



> 0, 



(49) 



where we identified uj with the bulk angular velocity. Eq. 49 is related to the conservation 



of energy. However, when condition 48 is not fulfilled, it gives rise to the magnetorotational 
instability (MIR)!™, that appears as the result of the interplay of its three different terms: 
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i) the angular momentum acquired by the fluid (or particles); ii) an interaction term due to 
the coupling between the fluid angular momentum with the driven angular velocity; iii) the 
fluid thermal energy and configurational entropy. 

A typical experiment consists in the rotation of a fluid between two concentric cylinders 
- related to the so called Taylor-Couette instability - driven by velocity gradients. In the 
presence of an axial magnetic field, the Taylor-Couette instability is onset when Eq. 49 is 



not satisfied. Also, we may expect that, owing to the fact that for two different species with 
different inertial moment, I a ^ it can be expected that at some point r = r c of the radial 



axis an inversion of the sign of the inequality of Eq. 48 must take place, and instability is 



onset. In particular, in the presence of two different species with different inertial moments, 
the fluid may well be intrinsically unstable at enough high angular speed. MIR instabilities 
concur against the stability of plasma configurations and in the 1960's, when MHD power 
plants were considered to be an economic process to convert thermal energy into electrical 
energy, E. Velikhov, in 1962, discovered the electrothermal instability which is at the source 
of strong magnetohydrodynamic turbulenc e! 35 * 37 * 38 ! 



The above Eq. 47 can be written in the form of a conservative equation for energy: 



0_ 
di 



.(a) • 



>) 



|JW|2 

2/( Q ) 



+ (w ■ J (a) ) - T6 (Q) 



|J(")| 2 
2/( Q ) 



+ (u ■ J (Q) ) - T& {a) 



(V r(Q) V Q )). 



(50) 



Finally, we can transform the above Eq. 50 into a Poynting's Theorem for rotating fluids: 



(51) 



after defining a kind of Poynting vector for rotational fluids, which gives the rate of rotational 
energy flow: 



s = - W° 

2 ^ 

a 



|J(«)|2 

2JW 



+ (u/ Q) • J (a) ) - T6 (q) 



(52) 



The power input driving the process (source/sink term) is given by: 



|JW| 2 
2J( Q ) 



-(u- 3 {a) )+TG ia) 



(53) 
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and the total energy is defined by summing up the different contributions: 



W = a + «P a +«P B + «p rot . (54) 

Here, the term represents the thermal energy associated to specie (a), equal to 

—AF, the free energy of the physical system. For a system in contact with a reservoir at 
constant temperature this is the maximum work which can be done, the free energy tends 
to decrease for a system in thermal contact with a heat reservoir. In particular, notice that 



when the angular velocity cu is multiplied by Eq. 40, it is obtained the driving power. It is 
worth to point out the presence of the term u ■ J which plays an analogue role to the slip in 
electrical induction motors, that is, the lag between the rotor speed and the magnetic field's 
speed, provided by the stator's rotating speed. Furthermore, we may notice that the power 
input depends on the fluid compressibility V • v. This means that compressibility is a factor 
that determines the amount of transported angular momentum through the stress-tensor 
Tij, and may be responsible of new driving mechanism in addition to the well-known MRI. 
The driving energy of the rotating system can be expressed in the form: 

Sdriv = ~ (w • J) - AF. (55) 
Zip 

Next, we will discuss several examples, illustrating the application of the variational 
method. 



A. Power extraction from an hurricane 

A hurricane is a natural airborne structure that converts all forms of kinetic and potential 
energy by means of the transport of angular momentum from the inner core to the outer 
regions, conveyed either directly by moving matter, or by non-material stresses such as those 
exerted by electric fields or magnetic fields 33 . We may apply Eq. p3^ to this specific problem 
and obtain: 

* , =«(£-- J -H = -ir (56) 



or 



V = oj£ driv . (57) 
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Let us consider specifically the case of a hurricane in an axisymmetric configuration, with 
J = ul. Likely we have u 2 I 2 /I p 3> u 2 I + AF. We can envisage a simple and rough model of 
a hurricane, with total mass M and radius R, turning like a solid cylinder, and / = MR 2 /2. 
Hence, the total power driving the hurricane is given by 

. lu 3 M 2 R i ,M 2 R 4 

p = 8T^r Kw Sr- (58) 

or, as a function of the fluid density p, in the form: 

-7T 2 F? 6 7~ 2 

V = ^ 3 ^P 2 (59) 



Our result shows the same kind of dependency as shown by Chow & Chey^, and, in par- 
ticular, it shows that the intrinsic inertial momentum of the particles constituting the fluid 
plays a substantial role. 



B. Radiactive heating of planetary atmosphere and net angular momentum 

It is experimentally showrP^ that periodic radiative heating of the Earth's atmosphere 
transmits angular momentum to it due to the Earth-atmosphere coupling through frictional 
interactions^. The images sent by ESA's Venus Express confirms this fact on Venus (Earth's 
twin planet), with the presence of a "double-eye" atmospheric vortex at the planet's south 
pole and high velocity winds whirling westwards around the planet completing a full rotation 
in four days. 

Schubert and Withehead's^ did an experiment with the purpose to provide an explana- 
tion for the high velocities of apparent cloud formation, in the upper atmosphere of Venus. 
In this experiment, a bunsen flame rotating below a cylindrical annulus, filled with liquid 
mercury, induced the rotation of the liquid mercury, in a counter-direction to that of the 
rotating flame. The speed of the flame was 1 mm/s, and the temperature of the mercury 
increased from room temperature, at the rate of about 3°C per minute. After 5 minutes, a 
steady-state flow was established, with the mercury rotating in the counter-direction of the 
flame, with a speed about 4 mm/s. Considering that the entire mercury liquid as a spinning 
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body, we can estimate from Eq. 49 that 



f 1 ^- - ^ + TG^j > => JA(u/ + - w) + AT6 > 0. 



dr 



(60) 



Hence, the following result is obtained: 



/A ^ 1AT6 . , 

w'Ao;>-- — . (61) 

Here, w is the angular speed imposed on the system; u' the induced angular speed, due 
to a sustained source of energy; Aw = uj' — u. We use the following tabulated data: 
& = 16.6 J.K -1 mol -1 for mercury at temperature of the experiment, p ~ 13.6 x 10~ 3 kg/m 3 , 
and consider that the volume of mercury is contained in the channel of the experimental 



apparatus, forming a rim with an average radius R = 30 cm. With Eq. 49 we estimate 
that, after one minute, the speed must be of the order 3.7 mm/s. It is clear that the sense 
of rotation, and its speed, depend on the latent heat stored in the planetary atmosphere 
and its temperature (through AT). Although the initial assumptions taken in the present 
formulation require further research, considering that interactive terms with the medium, 



conveyed for example through thermal diffusion coefficient, are not covered by Eq. 49, the 
agreement is reasonable, and offers a possible explanation for this effect. Once again, the 
problem of radiative atmosphere heating reveals the interesting interplay between energy 
and entropy, each one attempting a different equilibrium condition. 



V. CONCLUSION 



The information-theoretic method proposed in this paper constitute an alternate ap- 
proach applying the concept of maximizing entropy to the problem of out of equilibrium 
physical systems. It bears some resemblance with Hamiltonian formulation of dynamics, 
which expresses first-order constraints of the Hamiltonian H in a In dimensional phase 
space, revealing the same formal symplectic structure shared by classical mechanics and 
thermodynamics. Although the simplifying assumption of isothermal system rules out its 
ability to explain accurately such problems as, for example, the coherent transport of angular 
momentum in astrophysics, or a certain kind of laboratory devices (e.g., the Ranque-Hilsch 
effect), the present method open to understanding certain trends, in particular, predicting 
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theoretically for electromagnetic-gravitational vortex the angular momentum forced trans- 
port outwardly from a symmetry axis of rotation. The kind of Umov-Poynting theorem 
obtained is revealing of the interplay between entropy and energy, the energy tending to 
a minimum, and the entropy tending to a maximum, explaining the formation of physical 
structures. In particular, it was shown that compressibility is an important property in 
the transport of angular momentum, and a possible driving mechanism for instability. This 
development is seen to be advantageous and opens options for systematic improvements. 
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